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We briefly report our calculation of the 2-loop coefficient of the coupling constant renormalization 
function Z g in lattice perturbation theory. The quantity under study is defined through go = Z g g, 
where go (g) is the bare (renormalized) coupling constant. The 2-loop expression for Z g can be 
directly related to the 3-loop bare j3-function Pl(8o)- 

Our calculation is performed using overlap fermions and Wilson gluons, and the background 
field technique has been chosen for convenience. Our results depend explicitly on the number 
of fermion flavors (Nf) and colors (AO. Since the dependence of Z g on the overlap parameter p 
cannot be extracted analytically, we tabulate our results for different values of p in the allowed 
range (0 < p < 2), focusing on values which are being used most frequently in simulations. 
Plots of the 1- and 2-loop results for Z g versus p exhibit a non-trivial dependence on the overlap 
parameter. 

A longer write-up of this work may be found in Ref.[|ljj: 

http: //xxx.lanl . gov/abs/ 7 9 . 43 68 
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1. INTRODUCTION 

In later years, use of non-ultralocal actions which preserve chiral symmetry on the lattice has 
become more viable. The two actions which are being used most frequently are overlap fermions 
based on the Wilson fermion action and domain-wall fermions. Overlap fermions [^] are noto- 
riously difficult to study, both numerically and analytically. Many recent promising results from 
simulations with overlap fermions have appeared; see, e.g., Refs. Regarding analytical 

computations, the only ones performed thus far have been either up to 1 loop, such as Refs. [Jsj] — 
[11], or vacuum diagrams at higher loops [|12|, 13]. The present work is the first one involving 



non-vacuum diagrams beyond the I -loop level. 

We compute the 2-loop renormalization Z g of the bare lattice coupling constant go in the 
presence of overlap fermions. We relate go to the renormalized coupling constant g^ as defined 
in the MS scheme at a scale p. ; at large momenta, these quantities are related as follows 

a Ms(A) = Ok) + di{jxa)c^ + d 2 {}JLa)al + (1.1) 
(0Cq = go/471, a^jg = g^/4n, a : lattice spacing). The coefficient d\ has been known for a long 



time; several calculations of d 2 have also appeared, either in the absence of fermions |[14j, |15p, or 
using ultra-local fermionic actions [[1^, |l7| ]. Knowledge of d 2 (jXa), along with the 3-loop MS- 
renormalized j8 -function allows us to derive the 3-loop bare lattice j8 -function, which dictates the 
dependence of lattice spacing on go- Ongoing efforts to compute the running coupling from the 



lattice [ 18| , |19Q have relied on a mixture of perturbative and non-perturbative investigations. As an 
example, relating to asF (SF: Schrodinger Functional scheme), entails an intermediate passage 
through the bare lattice coupling and the conversion from a^j§- to CXq is carried out perturbatively. 

The paper is organized as follows: The theoretical background and the methodology of our 
perturbative lattice calculations are given in Section II. Section III describes the perturbative ex- 
pansion of the overlap action. Details on our computation, numerical results and plots of Z g can be 
found in Section IV. 



2. THEORETICAL BACKGROUND 



The renormalized j8 -function describes the dependence of the renormalized coupling constant 
g on the scale inherent in the renormalization scheme (chosen to be the MS scheme). A bare 
/3 -function is also defined for the lattice regularization (fiiigo)) 
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where a is the lattice spacing, ft the renormalization scale and gy^ (go) the renormalized (bare) 



coupling constant. In the asymptotic limit, one can write the expansion of Eq. (2.1) in powers of go 

Mso) = -bogl-b^gl-bUl-- (2-2) 
= - b oS 3 m -hg 5 m -b 2 g 7 m +... (2.3) 

The coefficients bo,b\ are well-known universal constants (regularization independent); bf (i > 2) 
(regularization dependent) must be calculated perturbatively. piigo) and Pigy^) can be related 
using the renormalization function Z g , that is 
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/ ,,dlnZ 2 \ 1 
/3 L (so)=fl-^^/j Z g P(g Q Z; 1 ) (2.4) 



Computing Z? to 2 loops 



Z 2 g (go,ap) = l+gl (26oln(aA) + z o) +£o (2*1 + + O(So) (2.5) 

and inserting it in Eq. ( |2.4| ), allows us to extract the 3-loop coefficient b\. All quantities in Eq. (|3|), 



except /i (for overlap fermions), are known. The expansion of Eq. (2A) in powers of g^ provides 
the following relation for b\ 

b^ = b 2 -b l lo + boh (2.6) 
Thus, the evaluation of b\ requires only a 2-loop lattice calculation, of the quantity l\ . 

The most convenient and economical way to proceed with the calculation of Z g (go,ap.) is to 

use the background field technique, in which the following relation is valid 

Z A (g ,a[i)Z*(go,ajl) = 1 (2.7) 

where Za is the renormalized function of the background field. In the lattice version of the back- 
ground field technique, the link variable takes the form 

Ufi (x) = e iagoQ » W • e iaA » w (2. 8) 

(Q^i- quantum field, A^: background field). In this framework, instead of calculating Z g , it suffices 
to compute Za- For the above calculation, we consider the background field 1PI 2-point function, 



both in the continuum and on the lattice; for notation and details, the reader can refer to Q14|]. The 
gauge parameter X must also be renormalized, in order to compare lattice and continuum results 

A=Z e Ao, Z Q = \+glz ( Q ] + ... (2.9) 
(Zq: renormalization function of the quantum field). Finally, Z g takes the form 



zl 



, , ( 2 - 10 ) 



1 +g 2 o (K 1 ' - v [l> ) +g 4 o «' - vW) + *o«8#^ 

(Ao = 1: bare Feynman gauge). Here, the only unknown quantity is v^ 2 \ the 2-loop 1PI lattice 
Green's function of the background field; its gluonic contribution (Nf=0) is known. One can find 
the expressions for the remaining, known quantities in Ref. [JXJ] _ 

The fermionic contributions are associated with the diagrams of Fig. 1 (v^ 1 )) and Fig. 2 
(y( 2 )). In the present work, v^ 2 ' is perturbatively calculated for the first time using overlap fermions 
and Wilson gluons. For completeness, we also compute the coefficient and compare it with 
previous results. For overlap fermions, the mass counterterm (the filled circle in Fig. 2) equals zero, 
by virtue of the exact chiral symmetry of the overlap action; consequently, diagrams 19 and 20 both 
vanish. Diagrams with infrared divergences become convergent only when grouped together (6+12, 
7+11, 8+18, 9+17). 





1 2 

Fig. 1: Fermion contributions to the 1-loop function v' 1 '. Dashed lines ending on a cross represent background gluons. 
Solid lines represent fermions. 
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Fig. 2: Fermion contributions to the 2-loop function v' 2 ' . Dashed lines represent gluonic fields; those ending on a cross 
stand for background gluons. Solid lines represent fermions. The filled circle is a 1-loop fermion mass counterterm. 



3. OVERLAP ACTION 

In recent years, overlap fermions are being used ever more extensively in numerical simula- 
tions, both in the quenched approximation and beyond. This fact, along with the desirable proper- 
ties of the overlap action, was our motivation to calculate the j3-function with this type of fermions. 
The important advantage of the overlap action is that it preserves chiral symmetry while avoiding 
fermion doubling. It is also 0(a) improved. The main drawback of this action is that it is necessar- 
ily non-ultralocal; as a consequence, both numerical simulations and perturbative calculations are 
extremely difficult and demanding (in terms of human, as well as computer time). 

The lattice overlap action is given by 

S merUp =a &1 £^(n)D N (n,m)^(m) (3.1) 

n,m 

where Djv(n,m) is the overlap-Dirac operator 



D N (n,m) = p 



S n ,, 



X- 



1 



X 



1 



D 



w - 



(3.2) 



Vxtx, 

and Dw is the Wilson-Dirac operator. The overlap parameter p is restricted by the condition < 
p < 2 to guarantee the correct pole structure of D^. The coupling constant is included in the link 
variables, present in the definition of X, and one must take the perturbative expansion of X in 
powers of go. This expansion in momentum space takes the form 

X(p',p) = Zo(p)(27i) 4 8 P (p' - p) +X l (p',p) +X 2 (p',p) +X 3 (p',p) +X 4 (p',p) +O(g 5 ) (3.3) 

V v ' * * ' V v ' 

tree— level I— loop 2— loop 

where Xo lS the inverse fermion propagator and X[ are the vertices of the Wilson fermion action (p 
(p f ): fermion (antifermion) momentum). Inside of X,'s there appear background and/or quantum 
gluon fields; their expressions are given in Ref. [[[]]. Upon substituting the expression for X,'s in 
the overlap vertices, the latter become extremely lengthy and complicated. For instance, the vertex 
with Q-Q-A-A-^- 1 ? consists of 724,120 terms. 

The perturbative expansion of leads to the propagator of zero mass fermions and to gluon- 
fermion-antifermion vertices (with up to 4 gluons for the needs of this calculation). After laborious 
analytical manipulations, the overlap-Dirac operator is expanded into terms with up to 4 gluons as 
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D N (k 1 ,k 2 )=D (k l )(2n) 4 8 4 (k l -k 2 ) + 'L(k l ,k 2 ) (3.4) 
Do(ki) is the inverse propagator and E(ki,k2) leads to the vertices with up to 4 gluons. Due to the 
large size of L(k\,k2), its explicit form is given in Ref. [|lj|. 

4. RESULTS 

For the algebra involving lattice quantities, we make use of our symbolic manipulation package 
in Mathematica, with the inclusion of the additional overlap vertices. The first step to evaluate the 
diagrams is the contraction among vertices; There follow simplifications of the color dependence, 
Dirac matrices and tensor structures. We use symmetries of the theory, or any other additional 
symmetry that may appear in particular diagrams, to keep the size of the expression down to a 
minimum. The external momentum p appears in arguments of trigonometric functions and the 
extraction of p dependence is divided into two parts: first we isolate terms that give single and 
double logarithms (a few thousand terms, expressible in terms of known, tabulated integrals) and 
then for the convergent terms we employ naive Taylor expansion with respect to p up to &{p 2 ). 
The coefficients of terms proportional to p 2 must be numerically integrated over internal momenta. 
The integration is performed by optimized Fortran programs which are generated by our Mathe- 
matica 'integrator' routine. Each integral is evaluated for several lattices with varying lattice size 
and for different values of the overlap parameter p . Finally, we extrapolate the results to infinite 
size lattices; this procedure introduces an inherent systematic error, which we can estimate quite 
accurately. Infrared divergent diagrams must be summed up before performing the extrapolation. 

We denote the contribution of the i th 1-loop Feynman diagram to v^\p) as V ( - (/?); similarly, 
contributions of 2-loop diagrams to v'" 2 \p) are indicated by V- 2 \p) 

,(0 , „2„2/ 7 ,(i) , t (2)lnflV\ , (41) 

(4.2) 



^{p)=N f \kPWp 2 {^+kf {4]t]2 



where p 2 = 4£^ sin 2 (/fy/2). The index i runs over diagrams, and the coefficients kf, c ; - depend 
on the overlap parameter p. Moreover, cf' = [c ( - ' l ' /N + c^'^N]. These coefficients must satisfy 
conditions coming from comparison with continuum results and usage of Ward Identities; we have 
checked that these conditions are verified by our results. 



£*F»=0, £cf)=0, £*f>4 Ecf) = I ^(3A.-I) (4.3) 



' l4i " cg> = — , cf) = V c$ = - 5 -N, cfj = ^l (4.4) 

15 3A f> 16 3 ' 17 3 > 18 3N 



In Table I we tabulate the total 1- and 2-loop coefficients = c^ 1 ' ^ = L/C- 1 ' 1 , 

C (U) = ^-c- 1 ' for several values of the overlap parameter (0< p <2), after extrapolating to L — > oo. 
The results for k^ 1 ' are in good agreement with corresponding results of Ref. [j8|]. In the Table and 
all Figures, the errors accompanying our results are entirely due to this extrapolation. In certain 
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cases with large systematic errors we extended the integration up to L = 46. In general, the overlap 
action leads to coefficients which are very small for most values of p . As a consequence, systematic 
errors, which are by and large rather small, tend to be significant fractions of the result for p > 1.4. 



p 


*W 


c (i,-i) 


c (i,D 


0.2 


0.01581702(2) 


-0.0044(1) 


0.0118(5) 


0.3 


0.0133504717(2) 


-0.00321(6) 


0.0045(1) 


0.4 


0.0116910952(1) 


-0.00244(4) 


0.0030(1) 


0.5 


0.0104621922(2) 


-0.00191(1) 


0.0022(6) 


0.6 


0.0095058191(2) 


-0.001606(6) 


0.00176(2) 


0.7 


0.00874441051(7) 


-0.001397(3) 


0.00145(1) 


0.8 


0.00813753230(4) 


-0.001241(1) 


0.00124(1) 


0.9 


0.00766516396(3) 


-0.001107(1) 


0.001051(9) 


1.0 


0.00732057894(3) 


-0.000979(1) 


0.000872(3) 


1.1 


0.00710750173(2) 


-0.000849(2) 


0.000710(8) 


1.2 


0.00703970232(7) 


-0.000706(3) 


0.00052(1) 


1.3 


0.0071425543(2) 


-0.000543(4) 


0.00033(3) 


1.4 


0.0074569183(2) 


-0.000335(7) 


0.00007(1) 


1.5 


0.0080467046(1) 


-0.00005(1) 


-0.0002(1) 


1.6 


0.0090134204(1) 


0.00034(1) 


-0.0004(1) 


1.7 


0.010526080(2) 


0.00093(6) 


-0.0021(5) 


1.8 


0.0128914(2) 


0.0020(1) 


-0.02(3) 



Table 1: Numerical results for the 1-loop (k^> = ^Jt\ ) and 2-loop (c 

i 

coefficients for different values of the overlap parameter p . 



(i-i) = W 1 - 1 ) ,.(1-1) = 



,(1,1) 



We plot the 1-loop coefficient k^ 1 ' (Fig. 3) and the 2-loop coefficients c^' ^, c^' 1 ' (Fig. 4) 
with respect to p. The extrapolation errors are visible for p < 0.4 and p > 1.7. 
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Fig. 3: Plot of the total 1-loop coefficient for different 
values of the overlap parameter p. 
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Fig. 4: Plot of the total 2-loop coefficients eC 1 ' 1 ), c ( - 1 ~ 1> > 
against p. 
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